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The electromagnetic response to microwaves in the mixed state of YBa2Cu30y(YBCO) was mea- 
sured in order to investigate the electronic state inside and outside the vortex core. The magnetic- 
field dependence of the complex surface impedance at low temperatures was in good agreement with 
a general vortex dynamics description assuming that the field-independent viscous damping force 
and the linear restoring force were acting on the vortices. In other words, both real and imaginary 
parts of the complex resistivity, p\, and p2, were linear in B. This is explained by theories for d-wave 
superconductors. Using analysis based on the Coffey-Clem description of the complex penetration 
depth, we estimated that the vortex viscosity 77 at 10 K was (4 ~ 5) X 10~^ Ns/m^ This value 
corresponds to lo^t ~ 0.3 - 0.5, where ljq and r are the minimal gap frequency and the quasiparticle 
lifetime in the vortex core, respectively. These results suggest that the vortex core in YBCO is not 
in the deeply superclean regime but in the moderately clean regime. Investigation of the moderately 
clean vortex core in high-temperature superconductors is significant because physically new effects 
may be expected due to d-wave characteristics and to the quantum nature of cuprate superconduc- 
tors. The behavior of Zs as a function of B across the first order transition (FOT) of the vortex 
lattice was also investigated. Unlike Bi2Sr2CaCu2 0j; (BSCCO), no distinct anomaly was observed 
around the FOT in YBCO. Our results suggest that the rapid increase of Xs due to the change of 
superfiuid density at the FOT would be observed only in highly anisotropic two-dimensional vortex 
systems like BSCCO. We discuss these results in terms of the difference of the interlayer coupling 
and the energy scale between the two materials. 



I. INTRODUCTION 

The properties of high-Tc superconductors (HTSC) in 
an applied magnetic field have revealed a rich variety of 
fascinating phenomena. In particular, investigations of 
the electronic structure in the mixed state of HTSC have 
been key experiments in recent condensed matter physics. 
In conventional superconductors (CSC) with s-wave gap 
symmetry, there are bound states of quasiparticles (QPs) 
localized in the vortex core. The level separation between 
the bound states, Ai?, is 



A£' EE huja - Ao/kpC, 



(1) 



where Aq is the bulk gap energy and kp and ^ are the 
Fermi wavenumber and the GL coherence length, respec- 
tively |]l|. Here, luq is regarded as the minimal gap fre- 
quency inside the core. 

In CSC, Tc is below 10 K and ^ - 100 A, which lead to 
AE < 0.1 K. This could not be resolved by available ex- 
perimental techniques. Indeed, in the scanning tunneling 
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spectroscopy (STM) experiments in CSC [|| ||] , the local 
density of states inside the core was shown as a zero en- 
ergy peak anomaly, which was interpreted as a smearing 
of the discrete levels in the core |^ . 

When describing the dynamics of QPs in the mixed 
state, the key parameter is the width of the level in the 
core, SE, which is connected with the scattering time, r, 
of the QPs in the vortex core. The relationship between 
SE and r was expressed as 



SE 



n 

T 



(2) 
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Since SE > AE in CSC, the vortices in CSC can be 
regarded as tubes with a normal metal of radius ^. In 
such a case, the vortex motion was well described with 
the Bardeen - Stephen (B-S) model 

In HTSC, such a simple description of the vortex core 
is incorrect for several reasons. First, the energy gap has 
d-wave symmetry in HTSC j?]. Since the amplitude of 
the gap at the node is zero, QPs are not localized in the 
vortex core but extended along the node direction ^, ||, 
p^ . A theoretical calculation suggested that there were 
no truly localized states in the vortex core in pure dx^-y^ 
superconductors On the other hand, recent STM 

results in HTSC showed that the peak at zero energy 
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vanished and that the density of states at a finite energy 
below the gap Aq was enhanced near the center of the 
vortex core [|l2|, |l3[ |4|. Therefore, the presence or the 
absence of locahzed states in the vortex core of HTSC is 
one of the central basic questions. 

The second different feature of the vortex core in HTSC 
is the semi-quantum nature of the core. Since Aq 20 
meV for HTSC, that means ^ - 20 A. This gives kp^ - 
20, which is quite different from the case of CSC (for 
which kpS, ^ 200, typically). Indeed, if the peak observed 
in the STM experiment mentioned above corresponds to 
the quantum levels in the core, it gives kpS. ~ 2. The 
physics of such a quantum core has been undeveloped up 
to now. 

Such anomalous aspects of the electronic structure 
around the core in HTSC should affect the dynamic re- 
sponse of the vortices. If the vortex motion is described 
phenomenologically, the dissipation can be discussed in 
terms of the frictional force. 



V^L + CtH^L X z. 



(3) 



which is proportional to the velocity of vortices, v^. Here 
T] and an are the viscous drag and the Hall coefhcients, 
respectively, and z is the unit vector parallel to an ap- 
plied magnetic field. The dissipation mechanism is, in 
turn, closely related to the characteristic scattering time 
of QPs in the vortex core, r. The vortex viscosity rj, and 
Hall constant an were calculated [|l5[ |l^ for arbitrary 
values of ujqt, and it was found that 



rj = irhn 
an — T^hn 



1 + ' 

1 + {loot)^'' 



(4) 
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where n is the QP concentration. The ratio of the level 
spacing AE to the level width SE, 



AE 

Je 



(6) 



is used as an index of the cleanness in the vortex core. 
In the dirty limit {loot < 1), the Hall effect (Eq. (|)) is 
negligible, and Eq. was expressed as 77 = nhnuJoT = 
Bc2^Q/Pm where p„ = m* /ne^r, and Bc2 is the upper 
critical field. This is the B-S expression of the viscosity, 
?7 1^. In the opposite limit, ojqt ^ l(the superclean 
limit), the Hall effect is dominant and rj is expressed as 
77 ~ Trhn/ujQT. According to the diagonal component 
of the vortex resistivity under constant transport current, 
J, was expressed as 



Pv 



where the effective viscosity r/cB was defined as 



(7) 



(8) 



This means that the experimentally obtained viscous co- 
efficient under the condition J = const., was r]cs. Equa- 
tions d), (|), and (|) yield 

?7cfr 



UJqT 



Trhn 



(9) 



Therefore, we can determine the QP lifetime in the vortex 
core from estimates of the viscosity of the vortex motion 
in the mixed state of HTSC. 

Due to the strong pinning in HTSC, dc-measurements 
are not appropriate to study the vortex dynamics at low 
temperatures. Instead, measurement of the microwave 
electromagnetic response is a powerful technique to probe 
the dynamics of the vortices in HTSC. From the high- 
frequency response of the vortices, several parameters 
relating to the vortex dynamics, such as the viscosity, 
pinning constant, and Hall constant can be estimated 
according to existing theories of vortex motion. 

There have already been many experimental ef- 
forts to determine the vortex dynamics parameters of 
YBa2Cu30j,, and a thorough review is provided in [pT| . 
To the best of our knowledge, experimental determina- 
tion of the vortex viscosity 77 in YBCO at low tempera- 
tures was first performed by Pambianchi et al., using the 
parallel plate resonator technique iQ. They measured 
the magnetic-field dependence of the complex resistivity 
up to 0.4 T at 11 GHz and showed that the tempera- 
ture dependence of 77 was consistent with a temperature- 
dependent normal state resistivity in the context of the 
B-S model A similar result was obtained in ||l^ by 
measuring the surface resistance Rs as a function of mag- 
netic field up to 8 T at several frequencies. A different 
report on the Rs measurement as a function of magnetic 
field up to 7 T at low temperatures [|o| , deduced a very 
large viscosity, which was two orders of magnitude larger 
than in conventional superconductors. This huge value 
was interpreted as evidence that the core of YBCO was 
in the superclean regime, ujqt 3> 1. However, their es- 
timation of the viscosity was based on the assumption 
that the pinning frequency LUp is much smaller than the 
measurement frequency uj. Under that assumption, they 
estimated the viscosity only from the data of the surface 
resistance, Rs (the real part of the surface impedance, 
Zs). In CSC, the assumption, uj ujp was reasonable 
in the microwave region since the characteristic pinning 
frequency ujp ~ 100 MHz [|l). In HTSC, however, the 
condition w ^ Wp might not be satisfied. Indeed, experi- 
mentally estimated values of the pinning frequency tUp at 
low temperatures in previous reports |jl^, ^ |2|] were on 
the order of 10 GHz. In such a case, the estimation of 
77 from Rs alone leads to an incorrect value for 77, and the 
measurement of Zs as a complex value is essential. Parks 
et al. |2^ measured the complex resistivity, p, of YBCO 
thin films as a function of frequency in the THz domain. 
They reported a relatively small 77 and insisted that the 
smallness was due to the effects of anisotropy and d-wave 
gap symmetry. If their interpretation is correct, there 
should be a definite difference in the magnetic-field de- 
pendence of the surface impedance from what would be 
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expected from the usual flux-flow theory. However, they 
reported only the frequency dependence of the complex 
resistivity p{uj) at a fixed field 6 T, and detailed field de- 
pendence data were unavailable. To discuss the deviation 
from the conventional vortex dynamics description, the 
magnetic-field dependence of p should be investigated. 

Considering the state of previous high-frequency mea- 
surements of 7y-values, there is a large ambiguity. To 
resolve such confusion, complex surface impedance mea- 
surements as functions of magnetic field, frequency, and 
temperature should be performed up to higher magnetic 
fields. 

In addition, it should be emphasized again that in 
HTSC, the vortex core is likely to be close to the quan- 
tum limit since kp^ is small and the major component 
of the condensate wave function has d-wave symmetry. 
In such a highly anomalous situation, to the best of our 
knowledge, there have been no rigorous calculations of 
the dynamic properties of vortices. So the appearance of 
some new effect may be expected. 

Another problem we focus on in this study is the rela- 
tionship between the electronic structure and the vortex 
structure in the mixed state. Previously, we found that 
in Bi2Sr2CaCu20j^, a rapid increase of surface reactance, 
Xs, took place at the first order transition (FOT) of the 
vortex lattice |2^. Since Xg is proportional to the real 
part of the effective penetration depth, it was proposed 
that this increase may be ascribed to the decrease in the 
superfluid density rather than to the loss of pinning. This 
suggested that the change in the structure of the vortex 
lattice largely affected the electronic structure of the vor- 
tex. If this phenomenon were observed in other cuprates, 
where the FOT occurred, it would be a universal new 
effect found in HTSC associated with the FOT of the 
vortex lattice. Thus, we have investigated how universal 
the phenomenon is in materials with different values of 
anisotropy. 

Furthermore, quite recently, the existence of a collec- 
tive mode in the microwave frequency range was pre- 
dicted in unconventional superconductors with mixed 
symmetry order parameters Such a mode might 

be expected to show up as a resonance in the surface 
impedance. 

To study the fascinating problems described above, in 
this paper, we present complex surface impedance mea- 
surements on untwinned YBCO single crystal as func- 
tions of temperature and magnetic field. These measure- 
ments have been made in the wide region of the mixed 
state (10 K < r < Tc, < < 15 T) at several fre- 
quencies, and provide information on electronic structure 
in the mixed state. In the following section, details of 
the experimental technique to measure Zs in an applied 
magnetic field together with the analysis based on vor- 
tex dynamics will be described. In Sec. Ill, the results of 
the Zg measurements will be shown. Both the real and 
imaginary parts of the complex resistivity show the B- 
linear behavior at low temperatures below 60 K. Through 
comparison with a theoretical calculation by Coffey and 



Clem |27| | , we estimate the pinning frequency, vortex vis- 
cosity and pinning constant. Estimated values of rj at 
10 K are 4 - 5 x 10^^ Ns/m'^. These values corre- 
spond to ljqt ~ 0.3 - 0.5, which suggests that the vortex 
core of YBCO is not in the deeply superclean regime but 
in the moderately clean regime. In the data from the 
high temperature region above 60 K, we cannot find any 
obvious anomaly around the FOT in both Rg and Xs- 
This is in contrast to the results observed in BSCCO. In 
Sec. IV, the observed i3-linear behavior is explained in 
terms of a theoretical consideration for d-wave supercon- 
ductors. The significance of the moderately clean vortex 
core realized in YBCO will be discussed together with 
the relationship between the Zs and the FOT in YBCO, 
and the difference between YBCO and BSCCO. More- 
over, we will briefly comment on the possibility of the 
existence of the collective mode predicted by Balatsky et 
ai. Sec. V is dedicated to our conclusions. 



II. EXPERIMENTAL 
A. Apparatus 

YBa2Cu30y single crystals were grown by a self- flux 
method using Y2O3 crucibles and were annealed at 
450°C. The sample described in this study was a nat- 
urally untwinned, slightly overdoped single crystal with 
dimensions of 0.6 x 0.6 x 0.05 mm^^. It showed a su- 
perconducting transition at Tc = 91.2 K, which was de- 
termined from the magnetization measurement using a 
commercial SQUID magnetometer (Quantum Design). 

The microwave surface impedance = i?s + iXs, 
where and X^ are the surface resistance and surface 
reactance, respectively, was measured by the cavity per- 
turbation technique with cylindrical cavities operated in 
the TEoii mode. Since a superconducting cavity with 
high sensitivity ||2^ could not be used in an applied mag- 
netic field, we used copper cavities tuned to 19.1, 31.7, 
and 40.8 GHz. The Q- values of each cavity at 10 K were 
- 6.5 x 10* at 19.1 GHz, ~ 2.8 x lO" at 31.7 GHz, and 

2.7 X 10^ at 40.8 GHz. The sample was located at the 
center of each cavity, at the antinode of the microwave 
magnetic field H,-{ being parallel to the c-axis of the sam- 
ple. Thus, the microwave currents flowed in the ab plane 
(Cu02 plane). Details will be described in a separate 
publication. Rg and Xg for the sample in zero field could 
be simultaneously obtained from the shift in the quality 
factor of the cavity, A(1/2Q), and the shift in the reso- 
nance frequency, A/. According to cavity perturbation 
theory ]2^ , in the skin depth regime where the microwave 
skin depth is significantly shorter than the dimensions of 
the sample, Rg and X^ may be expressed as 

X. = G(-A_A),c.G(-f)+C. ,11, 
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where Qs and fs are the quahty factor and the reso- 
nant frequency of the cavity including the sample, and 
Qb and /& are those of the empty cavity without the sam- 
ple. G is a geometrical factor that depends only on the 
shape of the sample and the cavity, and G is a metal- 
lic shift constant, which represents the shift of the res- 
onance frequency of the cavity with a perfect conductor 
with the same dimensions as the sample. In zero field, 
Xs is proportional to the London penetration depth, Al, 
as Xg = ploXl, where fi is the vacuum permeability and 
uj is the angular frequency. 

However, in Eq. (|l^), the effects of thermal expansion 
should be taken into account at high temperatures, as 
was pointed out in Therefore, Eq. ( [Tl| ) were cor- 

rected to the following expressions. 



Xs=G 



«(T) +G, 



(12) 



where a{T) is the temperature-dependent factor related 
to the relative difference of thermal expansivities between 
the cavity and the sample. In order to determine the ab- 
solute values of Rs and Xs , we should fix the two param- 
eters G, and G, and estimate a{T) in Eqs. ( p^ ) and (12). 
We determined G, G, and a{T), as follows. First, in the 
normal state, the relationship, Rs — (/icjp„/2)^/^ is ap- 
propriate, pn of the sample at 100 K is 6.4 x 10"^ m, 
which gives Rs = 0.28 ft at 31.7 GHz. Substitution of the 
value into Eq. ( p^ ) gives the value of G. Second, a{T) 
was estimated from the difference of the thermal expan- 
sivities between Cu [|ll, and YBCO |||. Finally, the 
value of G was determined by assuming that the London 
penetration depth at T = K, Al(0) ~ 1400 A Q. The 
corrected Xs data satisfied that (1) the Hagen-Rubens 
relation {Rs = Xs) was held in the normal state above 
Tc, and (2) at T ^ 0, Xs/puj was extrapolated to Al(0). 
After this correction, the temperature dependence of Xl 
in zero field was consistent with the data reported in . 
Rs and Xs in an applied magnetic field were subsequently 
determined by using the values of G, G, and a{T) in zero 
field. We confirmed that the magnetic field dependence of 
Zs was not affected by the corrections mentioned above. 
Surface impedance measurements were performed both 
with swept H after zero-field cooled conditions and with 
swept T under field-cooled conditions. In all measure- 
ments, the static applied magnetic field iJdc was parallel 
to the c-axis. 



B. Method of the Analysis 

In this subsection, we will explain the procedure used 
in this study to determine the vortex dynamics parame- 
ters. The complex surface impedance Zs is related to the 
complex penetration depth. A, as 

Zs — ipujX, (13) 

The complex resistivity p ~ pi + ip2 is expressed as 

p — Zs /ipoj — ipojX^ . (14) 



The general behavior of A in the mixed state of type- 
II superconductors was calculated by Coffey and Clem 
(C-C) ||2^. In this model, a rigid, massless vortex lattice 
with linear restoring force and viscous damping force was 
assumed. The basic equation of motion of vortex was 
described as 



riu. + KU — 00 J X 2 + i{t) 



(15) 



where u is the displacement of vortices, ry is the viscosity 
in the absence of flux creep, k is the pinning constant, J 
is the transport current density, 00 is flux quantum and 
z is the unit vector parallel to -ffdc- f(0 is the stochas- 
tic force, which corresponds to the effects of creep, and 
the Hall effect was neglected. Based on a self-consistent 
treatment of vortex dynamics including the influence of 
the normal component in the two-fluid sense, they de- 
rived an expression for the complex penetration depth 
A as a function of temperature, magnetic-field, and fre- 
quency. 



A^ 



1 (1 - e)r - s(e + r^). 



l + s2 (l + s2)(i+r2) 
s {e + r^) + s{l - e)r 



b . (16) 



Here, s = 2X\/5^j^ where Xl is the London pene- 
tration depth and o„y is the normal fluid skin depth 
— 2p„//icij). The normalized fleld 6, and normal- 
ized frequency r were expressed as, 

b = cjjoB/pujXlij - B/Br„ 
r = u/ujp, 



(17) 



(18) 



where iVp = {np/ri){h{v)Io{v)) / U/keT, 
U is an activation energy of the pinning potential, e = 
l//o(i^), and /„ is a modifled Bessel function of the first 
kind of the order n. 

At low temperatures, (T < 60 K), which corresponds 
to s — > and ^ oo, i.e. e ^ 0, effects of both the 
flux creep and the thermally activated normal fluid can 
be neglected, and the equation ( p^ ) becomes a rather 
simple expression. 



A^ 



1 + 



1 



(19) 



Rs and Xs normalized by (xluXl are deduced from 
Eq. (pj|), and expressed as 



Rs 



IiujXl 




b (2Q) 



5 




0.2 



b (21) 



Note that Eqs. ®, @, (|l|) contain only two pa- 
rameters, r and b. Variability of parameters r and b 
is limited within the range where C-C model is valid 
{Bel ^ B <C Bc2, where Bci and Bc2 are lower and 
upper critical field, respectively, and huj < Aq). There- 
fore, ranges of r and b are ^ r <C oo, and i3ci/-Bc2(~ 
0) ^ 6 <C Bc2/Bci{^ oo), respectively. 

Here we introduce "the impedance-plane plot" , where 
the horizontal axis is Rs/fJ-ujXL and the vertical axis is 
Xs/hloXl, as is shown in Fig. 0. Rg/fiujXL and Xg/fiLuXL 
of Eq. dg), and (||l|) were plotted through the parame- 
ter b with fixed r's. In the figure any curves in this plane 
show a concave upward curvature. Therefore, provided 
that the experimental data obeys the behavior in Fig |^, 
we can estimate the depinning frequency tUp uniquely by 
comparing the measured Zs{B) data with the curves in 
Fig. n^. A fter that we can determine 77 through Eqs. (20), 
and pl|). To sum up the procedure, first, we try to fit the 



data with the curve in Fig. ^ for a fixed r value. Substi- 
tuting the r value into Eqs. (^^ and (|l|) enables us to 
determine another parameter b. Finally, the comparison 
between b and experimental B gives the 77 value. 



III. EXPERIMENTAL RESULTS 

Figure ^ shows Zg as a function of magnetic field up 
to 15 T at several different temperatures. Here, the 
longitudinal axis of the lower panel in Fig. || (ARg) 
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FIG. 2: Magnetic fleld dependence of (a) Xs and (b) Ai?3 
of YBCO up to 15 T at 31.7 GHz at various temperatures. 
Continuous lines and open marks represent the results from 
swept magnetic-field measurements and swept temperature 
measurements, respectively. 



means the change from the zero-field value, defined as 
ARs = Rs{B) — i?s(0). Continuous lines and discrete 
open marks represent the results from swept H measure- 
ments and swept T measurements, respectively. Since 
both data sets were consistent with each other, extrinsic 
effects caused by an inhomogeneous field distribution in 
the data with swept H measurements did not affect the 
essential features of the data discussed below. 

In Fig. 1^, the magnetic field dependence of Rs and Xs 
is roughly proportional to B^^^ over the whole range, 
which is expected in the high field region of the C-C for- 
mula. To see this more clearly, we plot the data in terms 
of the complex resistivity, Api = pi{B) — pi(0), and 
Ap2 = P2{B) — p2(0). In Fig. ^, we plot the magnetic 
field dependence of Api and Ap2 at several fixed tem- 
peratures. If i?s(0) <^ Rs{B), Api, Ap2 were expressed 
as Api ~ 2ARsXs/puj, and Ap2 ^(A{Xf) - ARD/pu, 
respectively. According to Eqs. ( p^ ) and (p^, within 
the low temperature approximation, both Api and Ap2 
were proportional to a magnetic field, B. As shown in the 
lower two panels of Fig. g ((c), and (d)), both Api and 
Ap2 change linearly as a function of B at low tempera- 
tures below 60 K. These behaviors suggested that the C- 
C model was appropriate to describe the vortex dynamics 



in YBCO. However, at higher ternperatures above 70 K, 
as was clearly seen in Api (Fig. || (a)), there is a clear 
deviation from the i3-linear behavior. This deviation be- 
came more prominent with increasing temperature. In 
the high temperature region, thermal effects on the vor- 
tices cannot be neglected, and the low temperature ap- 
proximation will no longer accurate. To obtain precise 
values for the vortex parameters in this regime, analysis 
including flux creep, and/or the i?-dependent behavior 
of the vortex dynamics parameters is needed. 

In Fig. ^ (a), we plot the magnetic- field dependence of 
Zs shown in Fig. ^ in the impedance plane. Results from 
the measurements at 19.1 GHz up to 5 T, and those at 
40.8 GHz up to 7 T are also shown in Fig. ^ (b), and (c). 
All data show a clear concave upward behavior, which 
suggests that we can analyze the data in terms of the 
vortex dynamics. The fitted curve based on Eqs. ( |20| ) 
and (21), together with the resulting r value, are also 
shown in Fig. ^. Indeed, the result of the fit is very good 
in agreement. Thus, we can determine Wp and ?/, by the 
method described in the previous section. 

In Fig. |[ we show the temperature dependence of the 
viscosity, r] and the pinning frequency, cjp. r] decreases 
with increasing temperature. The observed temperature 
dependence of r\ was found to scale with that of the upper 
critical field, _Bc2- In fact, if we assume that the temper- 
ature dependence of Bcfi as Bc2 — Bc2(0)(l — i^)/(l 
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FIG. 3: Magnetic field dependence of (a) Api at higher tem- 
peratures, (b) Ap2 at higher temperatures, (c) Api at lower 
temperatures, and (d) Ap2 at lower temperatures, at 31.7 
GHz. The origin of each curve has been shifted arbitrarily. 
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FIG. 4: (a) The impedance plane plot of the data at 31.7 
GHz at various temperatures. The solid lines represent fitting 
curves to the vortex dynamics model (Eqs. (^^ and (pl|)), 
together with the resulting r values. Actual magnetic-field 
values of each line are shown as dotted lines (5, 10, and 15 
T, respectively). Similar plots based on the data at 19.1 GHz 
up to 5 T, and at 40.8 GHz up to 7 T are shown in (b), and 
(c), respectively. 



where t = T/T^, and Bc2(0) = 80 T and p„ - 3 x 
10~^ ilm, the solid curve in Fig. |^ (a) was obtained. 
The data are almost in agreement with this approximate 
curve. Deviations might be caused by the temperature 
dependence of the normal state resistivity, p„, and/or 
an ambiguity in the high temperature region that will 
be discussed later. The pinning frequency ujp decreased 
with increasing temperature due to thermal effects. Such 
behavior in the temperature dependence of ?? and u]„is al- 
most consistent with the previous reports ||l8| , p^ , p3| . 
As was already mentioned, Eqs. (^0|) and ( |2lD are approx- 
imate expressions in the low temperature region, and in 
fact, the fit was not so good above 70 K. Therefore, the 
absolute values of these parameters at high temperature 
(above 70 K) might have some ambiguity. However, in 
the low temperature region, on which we mainly focus in 
this study, these parameters can be well defined within 
this approximation. 

It was surprising that the obtained values of rj and tOp 
from the measurements at different frequencies are al- 
most the same. This fact indicates that the frequency 
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at the FOT in YBCO within our experimental resolu- 
tion. This result is in contrast to what was observed in 
BSCCO with the most remarkable being the anisotropy. 
One possible reason for the absence is that the change 
of London penetration depth, Al, is too small to be de- 
tected. In YBCO, the anisotropy parameter 7 is smaller 
than that of BSCCO and the phase transition occurs at 
rather higher fields, which are the order of several tesla. 
In such high fields, the density of vortices is so large that 
the contribution of moving vortices to Zs may hide the 
subtle change of the superfluid density. However, if the 
anomaly is supposed to be as the same magnitude as in 
BSCCO, it should be AX^/X^ - 4 - 5% , which is weU 
above the measurement noise level (< Imil, ^ 1%). We 
will discuss the reasons for this difference later. 



IV. DISCUSSION 

We found that the change in the resistivity was linear 
in B. We also clarified that the flux-flow in YBCO was 
well described in terms of the Coffey-Clcm unified the- 
ory of vortex motion with _B-independent, and frequency- 
independent vortex dynamics parameters, and essentially 
not different from that in CSC. As will be seen below, this 



FIG. 5: Temperature dependence of the viscosity, r; (a), and 
the pinning frequency, ujp (b), estimated from the analysis 
(see text). Solid curve in the panel (a) are a approximate 
expression, r?(r) = ??(0)(1 - + t^). 



dependence of Zg was also in good agreement with the 
vortex dynamics description with frequency-independent 
values for 77 and ujp. The viscosity rj and the depinning 
frequency ujp at 10 K were found tobe~4-5 x 10^^ 
Ns/m^, and ~ 40 - 50 GHz, respectively. Since we kept 
J = const, in our experiments (where J is the amplitude 
of the microwave current), the viscosity obtained from 
the data corresponds to rjcs in Eq. (H). In YBCO, the 
carrier concentration n was proposed to be n = 0.25 per 
Cu02 plane |Q. By using Eq. (||), we estimate that ujot 
is ~ 0.3 - 0.5. Since the ujp was found to be 40 - 50 GHz, 
even in the measurements in the microwave and millime- 
terwave regions, we cannot conclude that thepinning is 
negligible. In that sense, the analysis in Ref. [Q, leading 
to the conclusion of the superclean core, was incorrect. 
Our results confirmed that the measurement of both the 
real part and the imaginary part is essential for HTSC. 

Next, let us move to the behavior of Zg around the 
FOT in the vortex lattice. Figure ^ shows the magnetic- 
field dependence of Zg, Api, and Ap2 up to 15 T at 80 
K, and 85 K. The dashed lines in each panel show the 
field corresponding to the FOT, which was determined 
by the resistivity data of a sample with the same dop- 
ing level p6|. We could not find any obvious anomaly 




5 10 150 5 10 15 

Applied Field (T) Applied Field (T) 



FIG. 6: Magnetic field dependence of Zs, Api, and Ap2 
around the FOT of the vortex lattice. Left three panels ((a), 
(b), and (c)) and right three panels ((d), (e), and (f)) are 
the data at 80 K, and 85 K, respectively. Dotted lines in each 
panel represent the FOT field at each temperature, which was 
estimated from the resistivity data of sample from the same 
batch. 
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is consistent with theoretical considerations for d-wave 
superconductors ^ 

Theoreticahy, flux-flow was first discussed by Bardeen 
and Stephen. Q . They calculated the flux flow resistivity 
in the dirty limit, SE ~ ?i/r > Ag for s-wave supercon- 
ductors. After that, Larkin and Ovchinikov |^ showed 
results for the case of a moderately clean core with s-wave 
gap symmetry, and Kopnin and Kravtsov [ p^ calculated 
the conductivity and the Hall angle taking the Hall effects 
into account. Recently, Kopnin and Volovik |37| , and 
Makhlin |Q described the flux flow in the clean, d-wave 
case. They showed that, in the moderately clean limit, 
the QPs outside the core did not make an important con- 
tribution to the flux-flow conductivity. The expression 
for the flux-flow conductivity in a d-wave superconduc- 
tor was similar to that for an s-wave superconductor in 
the moderately clean regime. This means that the dis- 
sipative process mainly occurs at the boundary of the 
core, and is independent of the details of the electronic 
structure inside the core in the clean regime. Again, our 
results suggest that the flux-flow resistivity of YBCO was 
_B-linear, which is not so much different from that in a 
conventional s-wave superconductor. This result is quite 
consistent with the theoretical description of the flux-flow 
in a d-wave superconductor. 

Next, we found that the vortex viscosity r] was 4-5 
X 10"'' Ns/m^ at 10 K, which means that the vortex 
core of YBCO is in the moderately clean regime. As was 
described earlier, the electronic structure of vortices in 
HTSC is unusual because of (1) the d-wave symmetry of 
the energy gap, and (2) the quantum nature of the lev- 
els in the core. However, most of the existing theories 
took the quasiclassical approach to calculate the elec- 
tronic structure of the core. In HTSC, as kp^ is ~ 10 
and hcoo is ~ 10 K, the vortex core is likely to approach 
the quantum limit {kp^ ^1). In the quantum limit, since 
the minimal gap hu/Q is nearly equal to the energy gap Aq, 
the distinction between QPs inside and outside the core 
becomes meaningless. Larkin and Ovchinikov ]4C| |, and 
Koulakov and Larkin calculated the I-V character- 
istics and microwave absorption in the quantum, clean, 
s- wave cases, respectively. However, to our knowledge, 
there was no theoretical prediction on the flux flow state 
in the quantum, moderately clean, d- wave case. All of 
these estimations mean that the vortex core in HTSC is 
in the marginal region from the different points of view 
discussed above. 

Next, we will discuss the absence of any anomaly in the 
microwave response in YBCO at the FOT. In BSCCO, 
the rapid increase of Xs at the FOT was observed in 
both optimally- and over-doped samples Q. Note that 
no obvious anomaly was observed in Rs at the FOT. This 
provides strong evidence that the origin of the anomaly 
could be considered as arising from the change of the 
superfluid density. 

One possibility for the absence of the anomaly is that 
the mechanism of the FOT is different between BSCCO 
and YBCO. A key difference between YBCO and BSCCO 



is anisotropy. According to |Q, the characteristic field, 
Bcr, which separates the regions of three-dimensional 
and two-dimensional vortex behavior, was expressed as 
Bcr ^ 'Po/l'^'Pj where 7^ = rric/mab {iTiab and rric are 
the anisotropic effective mass in the ab plane and along 
the c axis, respectively), is the anisotropy parameter and 
d is the distance between adjacent Cu02 layers. Typical 
values for 7^ are 50 and ~ 20,000 in YBCO and BSCCO, 
respectively. While the vortices in BSCCO have a 2D-like 
character, the vortex lattice in YBCO is rather 3D-like. 
From this point of view, our results suggest that the rapid 
increase of Xs due to the change of superfluid density at 
FOT would be observed only in a highly anisotropic two- 
dimensional vortex system like BSCCO. In BSCCO, it 
has been established that the interlayer coupling changes 
dramatically at the FOT |44 . This change of interlayer 



coupling can cause a change in the current distribution 
in the sample on a macroscopic scale, which possibly pro- 
duces the change in the in-plane properties such as Rs 
and Xs However, if it is the case for the anomaly 
we have discussed, a similar anomaly should appear also 
in Rs, which is quite different from what we observed. 
Thus, we do not believe that such a change is responsible 
for the anomaly. 

Another possibility is the difference of energy between 
the two materials. In YBCO, the complex resistivity as 
a function of magnetic field was measured through the 
FOT fine up to 50 MHz |46). They reported definite 
changes in both real and imaginary parts of the com- 
plex resistivity, and interpreted the results as an abrupt 
change of the viscosity and the collapse of the shear mod- 
ulus due to the melting of the vortex lattice. However, 
in the microwave response, no sign of the FOT was ob- 
served. We speculate that the signal due to the change 
of the electronic state, may exist in another frequency 
window, lower or higher than the microwave range. 

Finally, we briefiy mention the attempt to detect the 
new resonance. Quite recently, Balatsky et al. predicted 
the existence of a new collective mode ||2^ whose fre- 
quency is tunable by magnetic field. According to [ p6[ , 
LUci ^ 72 X H/Hc2 GHz. Assuming that Hc2 ~ 100 T, 
LUci at 15 T is ~ 11 GHz, which is lower than our mea- 
surement frequencies. If Hc2 shows a usual BCS-like T 
dependence, the resonance originating from this mode 
would exist in a rather high T region below 15 T. How- 
ever, we could not observe such as a resonance structure 
in our measurements of Zs{B), and Zs{T). If the life- 
time of this mode were strongly influenced by thermal 
smearing, this mode would be observed in the rather low 
temperature region. 

In order to investigate the behavior around the FOT 
and find the possible collective excitation in HTSC, a 
measurement at several GHz under high magnetic field 
is now in progress. 
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V. CONCLUSION 

The microwave electromagnetic response in the mixed 
state of YBa2Cu30j, (YBCO) was measured in order to 
investigate the electronic state around the vortex core. 
The magnetic-field dependence of the complex surface 
impedance at low temperature was in good agreement 
with a general vortex dynamics description assuming that 
the field-independent viscous damping force and the lin- 
ear restoring force acted on vortices. In other words, both 
the real and imaginary parts of the complex resistivity, 
pi, and p2 were linear in B. This is explained by the- 
ories for d-wave superconductors. From analysis based 
on the Coffey-Clem description of the complex penetra- 
tion depth, we estimated that the vortex viscosity 77 at 
10 K was 4 - 5 X 10^'' Ns/m^. This value corresponds 
to wqt ~ 0.3 - 0.5, where wq and r are the minimal 
gap frequency and the quasiparticle lifetime in the vor- 
tex core, respectively. These results definitely show that 
the core of YBCO is in the moderately clean regime. 
Investigation of the moderately clean core of HTSC is 
significant because physically new effects are expected 
from the characteristic rf-wave and quantum nature of 



cuprate superconductors. The behavior of Zs as a func- 
tion of B across the FOT was also investigated. Unlike 
in Bi2Sr2CaCu20y (BSCCO), a distinct anomaly was not 
observed around the FOT in YBCO. Our results suggest 
that the rapid increase of Xg due to the change of su- 
perfluid density at the FOT would be observed only in 
a highly anisotropic two-dimensional vortex system like 
BSCCO. 
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